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ABSTRACT  (CmNmm  on  nfwii  N*  II  !#• muff  Sp  Mack  m—Sarj 

ASL-DUST  Is  a  computer  code  that  models:  (1)  the  dust  clouds  generated  by 
tactical  munitions  Impacting  the  ground,  and  (2)  the  propagation  of  electro- 
optical  signals  through  the  resulting  dusty  environment.  This  volume  discusses 
the  models  In  ASL-DUST  and  compares  the  code  results  with  the  experimental  data 
base.  Volume  2  is  the  User's  Manual. 


PREFACE 

Several  efforts  to  model  the  general  problem  of  battlefield  obscuration, 
and  specifically  that  of  battlefield  dust,  have  been  carried  out  at  and 
under  the  auspices  of  the  Atmospheric  Sciences  Laboratory  (ASL)  over  the 
past  two  years.  One  such  effort  Is  documented  here;  namely,  the 
development  of  a  general  research  conputer  code  to  simulate  the  effects 
of  artillery  produced  dust  on  the  transmission  of  electromagnetic 
energy.  The  purpose  of  this  code  Is  to  assist  the  modeler  In  the 
Interpretation  of  past  test  data,  and  to  aid  In  the  methodology  and 
structure  of  future  tests.  As  Is  always  the  case  with  an  on-going 
effort,  this  report  represents  a  snapshot  In  time  of  a  dynamically 
evolving  code.  It  Is  necessary  to  modify  the  code  as  the  year 
progresses,  and  elements  of  this  effort  will  finally  be  Incorporated 
Into  the  Electro-Optics  Systems  Atmospheric  Effects  Library  (EO-SAEL). 


Melvin  6.  Heaps 

US  Army  Atmospheric  Sciences  Laboratory 
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SECTION  1 
INTRODUCTION 


A  wide  variety  of  electro-optical  (E-O)  sensors  are  employed  in 
the  modern  tactical  battlefield  environment.  Dust  clouds  are  generated 
when  tactical  munitions  detonate  at  or  below  the  ground  surface.  These 
dust  clouds  can  be  a  major  source  of  degradation  for  the  battlefield 
performance  of  the  E-0  sensors. 

In  Reference  1,  our  previous  report,  we  began  the  task  of  modeling 
the  munition  dust  clouds  and  the  propagation  of  E-0  signals  through  the 
clouds.  Under  the  present  contract  we  have 

•  Continued  the  dust  cloud  model  development 

•  Written  a  computer  program,  ASL-DUST,  incorporating  the 
models  developed  for  the  dust  clouds  and  the  propagation 
of  E-0  signals  through  the  clouds 

•  Validated  the  models  by  comparing  model  predictions  with 
experimental  test  data. 

This  volume  presents  the  new  dust  cloud  model  developments  and  the 
comparisons  with  test  data.  Volume  2  (Reference  2)  is  the  User's 
Manual  for  the  ASL-DUST  computer  code. 
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SECTION  2 

NEW  AND  IMPROVED  DUST  CLOUD  MODELS 


In  this  section  we  continue  the  modeling  effort  begun  in  Reference- 
1.  The  new  and  improved  dust  cloud  models  include 

•  Fractionization,  that  is,  mixing  of  particles  among  the 
different  size  groups 

•  A  base  cloud  model;  this  is  the  ground  level,  nonrising 
base  surge  dust  region 

•  Time  delays  before  horizontal  transport  of  the  main  and 
base  clouds  due  to  wind  transport 

•  A  vertical  profile  for  the  mean  wind  speed 

•  An  improved  fallout  model  for  the  main  dust  cloud. 

FRACTIONIZATION 

We  divide  the  particle  size  distribution  into  a  number  of  size 
groups.  The  rise,  wind  transpor* ,  diffusion,  and  fallout  of  each  size 
group  are  calculated  separately.  In  the  preliminary  model  we  assumed 
that  all  particles  of  a  given  size  group  remained  within  their  volume; 
we  assumed  nonmixing  of  particles  from  one  size  group  volume  to  another. 
At  later  times,  differential  fallout  would  produce  an  absolute  separa¬ 
tion  in  the  cloud,  with  the  lightest  particles  at  the  top  of  the  cloud 
and  the  heaviest  particles  at  the  bottom.  In  a  real  dust  cloud  we 
would  expect  to  find,  on  the  average,  lighter  particles  at  the  top  and 
heavier  particles  at  the  bottom,  but  not  the  perfect  separnMo  of  the 
preliminary  model. 

There  are  many  processes  in  the  rising  duct  cloud  that  tei  .0 
mix  the  different  size  particles,  especially  the  lighter  particles 


among  the  heavier  particles.  Turbulence,  the  circulating  air  flow 
fields  of  the  rising  cloud,  nonuniformity  of  the  flow  fields  over  the 
volume  of  the  size  group,  and  adhesion  of  smaller  particles  to  each 
other  and  to  the  larger  particles  will  all  produce  mixing.  Since 
heavier  particles  are  less  affected  by  the  mixing  processes  than 
lighter  particles,  we  expect  mixing  to  be  more  efficient  as  the 
particle  size  decreases. 

We  formulate  the  following  simple  model  to  simulate  the  particle 
fractionization.  We  assume  that  a  fraction  of  the  particles  of  a 
given  size  group  is  mixed  among  all  larger  size  groups.  Thur  a 
fractionized  size  group  will  consist  of  a  mixture  of  particles  from 
itself  and  from  all  smaller  size  groups.  We  take  a  constant  fraction 
of  the  mass  from  each  size  group  and  distribute  this  mass  among  all 
larger  size  groups  in  proportion  to  the  mass  fractions  of  the  heavier 
groups. 

Let 

Fz  =  fraction  of  the  mass  of  a  size  group  that  is  mixed  among 
the  heavier  size  groups 

F^  =  mass  fraction  of  the  unmixed  size  group  i  (mass  of 
particles  in  size  group  i  before  mixing  divided  by 
total  mass  of  all  particles) 

pu  =  propagation  mass  coefficient  (for  extinction,  absorp¬ 
tion,  scatter,  or  backscatter)  for  unn.i.xed  size  group 

.  ,  2  -1, 

1  (cm  g  ) 

n 

F,^  =  ^  =  sum  of  the  mass  fractions  of  all  size  groups 

j  =  i+l  larger  than  i.  Note  that  F,^  is  only  defined 

for  i  between  1  and  n-1,  where  n  is  the  total 
number  of  size  groups. 


The  mass  contributed  to  size  group  j  (j  >  i)  from  size  group  i  after 
mixing  is 
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where  M  is  the  total  mass  of  all  particles  (g) .  The  total  mass  in 
size  group  j  after  mixing  is  the  mass  of  size  group  j  that  is  not 
contributed  to  larger  groups  plus  the  sum  of  all  mass  contributions 
to  group  j  from  size  groups  smaller  than  j: 


'MFMia-FJ 


j  =  1 


i-1  F. 


MF...S.  +  MF  . F  > 

•  Mj  i  Mi  z  f-1. 


2  <  j  <  n 


i=l  Ti 


where 


M.  =  mass  in  size  group  j  after  mixing  (g) 
(l  -  vz  i<n 


VI  j  =  n 

The  factor  represents  the  fraction  of  unmixed  size  group  j  left 
after  contributing  mass  to  larger  groups.  The  last  group  (j  = n)  has 
no  larger  groups  and  thus  does  not  contribute  mass. 

After  mixing,  each  size  group  has  a  different  mass  and  (except 
for  size  group  1)  a  different  group  propagation  mass  coefficient.  The 
mass  , raction  of  size  group  j  after  mixing  is 
M.' 

FMj  "M*  (4) 

The  group  propagation  mass  coefficient  after  mixing  is 


j  =  1 


L1  F. 


2  -1 
cm  g 


FM.S.p.  +  FU.F  £  y. 

Mj  J  J  MJ  z  FT.  Mi 


2  <  j  <  n 


We  assume  that  the  rise,  transport,  and  diffusion  of  each  si  c 
group  is  still  given  by  the  models  developed  for  the  unmixed  size 
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groups.  Hence  the  cloud  location  and  shape  are  not  changed;  only 
the  population  of  particles  within  the  cloud  is  mixed.  The  amount  of 
mixing  is  controlled  by  the  mixing  fraction  ?z .  For  the  present  we 
take 


A  value  of  one-half  allows  a  moderate  amount  of  mixing.  With  this 
value  model  predictions  of  electro-optical  transmissions  through  dust 
clouds  are  in  good  agreement  with  measured  test  results. 

BASE  CLOUD 

Figure  1  is  an  artist's  sketch  of  the  dust  regions  produced  by  a 
munition  impacting  the  ground.  In  Reference-  1  we  developed  models  for 
the  main  cloud  dust  region.  In  this  report  we  will  develop  models  for 
the  base  cloud  dust  region.  That  part  of  the  stem  due  to  fallout  from 
the  main  cloud  is  included  in  the  main  cloud  models.  But  no  separate 
geometric  stem  region  is  currently  modeled. 

When  a  munition  detonates  at  or  below  the  ground  surface,  a  quan¬ 
tity  of  soil  is  ejected  into  the  air  from  the  resultant  crater.  A 
region  of  rising,  circulating  air  flow  is  set  up  by  the  shock  wave  and 
the  buoyant  rise  of  the  heated  air.  Part  of  the  ejected  soil  is  en¬ 
trained  into  the  rising  flow  fields  and  forms  the  main  dust  cloud. 

Part  of  the  soil  ejecta  is  thrown  to  the  side  and  is  not  entrained. 
Also,  the  explosion-produced  shock  will  scour  some  additional  dust 
from  the  ground  surface.  The  base  cloud  dust  region  represents  that 
portion  of  the  lofted  dust  whi  h  is  nonrising. 

We  assume  that  the  initial  base  cloud  forms  instantaneously  at 
burst  time,  like  the  main  cloud.  We  use  the  same  spheroidal  geometry 
and  Gaussian  mass  distribution  as  in  the  main  cloud  model,  and  we 
scale  the  initial  radii  of  the  base  cloud  in  terms  of  the  initial 
radii  of  the  main  cloud.  Let 

R^  =  initial  horizontal  radius  of  the  main  cloud  along  the 
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Figure  1.  Munition-produced  dust  regions. 


artillery  shell  track  direction  (ml 
Rp  ~  initial  horizontal  radi  <s  of  the  mu  .n  cloud  in  the 
direction  perpendicular  to  the  shell  track  (m) 

Ry  -  initial  vertical  radius  of  the  main  cloud  (m) . 

We  expect  the  initial  vertical  dimension  of  the  base  cloud  to  be 
about  the  same  as  the  main  cloud,  but  expect  the  dimensions  to  the 
side  to  be  larger  due  to  the  ejecta  thrown  sideways.  We  take 


R 


BT  =  BMRT 


(6a) 


RBP  =  bmrp 


rbv  =  Rv 


(6b) 

(6c) 


where  RgT,  Rgp,  Rgv  are  the  corresponding  base  cloud  radii  in  the 
shell  track,  cross  track,  and  vertical  directions.  Bu  is  a  constant 


multiplying  factor.  Considering  the  base  cloud  formation  mechanisms, 
we  would  expect  to  be  in  the  range  of  about  2  to  4.  We  obtain  the 
best  agreement  between  model  and  test  results  by  choosing 


BM  =  3* 

We  assume  the  base  cloud  has  the  same  dust  and  carbon  composition 
as  the  main  cloud,  and  thus  we  use  the  same  size  groups.  At  burst 
time  all  size  groups  occupy  the  same  volume.  After  burst  time  the 
groups  begin  falling  out  of  the  cloud,  and  the  wind  transports  them 
horizontally.  Let  the  wind  be  along  the  positive  (horizontal)  X- 
axis,  Z  is  positive  upward,  and  Y  is  in  the  cross-wind  direction. 

Then  the  coordinates  of  the  centroid  of  the  smallest,  essentially 
massless  particles  are 


where 

V 


Xco<l>  =  VW  At 


Y  (t) 
co  J 


CO 


0 

(t)  =  R 


(7a) 

(7b) 


BV 


W 
At 


=  mean  wind  velocity  at  altitude  Zcq  (m  s-1) 

=  length  af  time  the  wind  has  been  transporting  the  base 
cloud  horizontally  (s). 


Since  the  base  cloud  is  nonrising,  the  Z  coordinate  of  the  massless 
particles  remains  at  the  initial  position. 

For  the  centroid  coordinates  of  the  finite  mass  size  groups,  we 
use  the  same  fallout  and  transport  models  as  the  main  cloud: 


;i(t)  =  Xco(t)  -±[in{(l.^) 


C2At 


•] 


Y  .  (t)  =  0 
ci v  ' 


(8b) 


Z  .  (t) 
civ  1 


Z  (t)  -  V. .At 
co^  '  ti 


(8c) 
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where 


X^(t)  =  X  coordinate  of  the  centroid  of  the  ith  size  group  (m) 

V  .  =  terminal  fall  velocity  under  gravity  of  the  ith  size 

Z1  -1 
group  (ms  ) 

and  Cj  and  are  the  same  main  cloud  constants  as  before. 

However,  the  diffusion  of  the  base  cloud  differs  from  that  of  the 
main  cloud.  During  the  rise  phase  the  main  cloud  diffusion  is  driven 
by  the  buoyant  rise  velocity  and  only  later  is  controlled  by  atmospheric 
diffusion  processes.  Since  the  base  cloud  is  nonrising,  its  diffusion 
is  atmospherically  controlled  from  the  beginning.  Table  1,  with  data 
taken  from  Reference  3,  shows  the  scaling  factors  for  angular  lateral 
spread  and  vertical  spread  due  to  atmospheric  diffusion  as  a  function 
of  the  atmospheric  Pasquill  stability  category.  The  relations  between 
the  scaling  factors  and  the  particulate  Gaussian  standard  deviations 
are 

4 . 3oy (km) 

0frad)  -  -Xiao-  (9) 

(^(km)  =  aX(km)S  ,  (10) 

where 
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TABLE  1.  ATMOSPHERIC  DIFFUSION  SCALING  FACTORS 


a,s  =  scaling  factors. 

Rewriting  the  scaling  formulas  in  terms  of  radii  (r = 2.15o)  and 
converting  distance  units  to  meters,  we  have 


9 (deg)  Vwt 

ry(t)  "  2(57.296)  m 

(11) 

rz(t)  =  2.15(a)103(1_s)[Vwt]S  m  . 

(12) 

These  radii  relations  are  for  a  point  source.  For  our  case  of  instan¬ 
taneous  finite  source,  we  take  the  diffusion  relations  for  the  smallest 
particles  to  be 

Rxo(t)  =  Rx(0)  +  ry(t)  m  (13) 

Ryo(t)  =  Ry(o)  +  ry(t)  (14) 

Rzo(t)  =  (Rz/S(0)  +  rz/S(t))S  ’  (15) 

where  Rx(o),  Ry(o),  and  Rz(o)  are  the  initial  base  cloud  radii  in 
the  x,  y,  and  z  directions,  respectively. 

For  the  finite  mass  size  groups,  we  assume  the  same  relation  as 
for  the  main  cloud,  namely 
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6V  .  \] 

#*') 


where  Rq  is  the  radius  of  the  zero  diameter  particle  group  and  is 
the  radius  of  the  ith  size  group.  Differentiating  and  rearranging, 
the  equations  for  the  radii  for  size  group  i  become 


R  .  (t)  = 
xi  J 


R  (t) 
xo  ' 


H 

t.T  ' 


R  (t)  = 
yi  1 


R  (t) 
yol  ' 


/ 16V  .  \ 

(-#•■) 


R  .  (t)  = 
ziv 


R  (t) 
:.o  J 


fH 


s/2  (2s  -  J  ) 


The  total  dust  and  carbon  mass  lofted  into  the  base  cloud  is  only 
a  small  fraction  of  that  in  the  main  cloud.  Typical  estimates  of  the 
base  cloud  mass  are  about  5  to  10  percent  of  the  main  cloud  mass. 

Many  electro-optical  sensor  sight  paths  in  the  battlefield  are  within 
the  first  2  or  3  meters  above  the  ground.  The  nonrising  base  cloud 
can  produce  significant  attenuation  effects  on  these  low-altitude 
sight  paths.  The  lighter  particles  in  the  main  cloud  rapid]  •  rise 
above  these  low  sight  paths.  Thus  even  though  the  main  clout  h«s  a 
much  larger  mass,  at  later  times  the  attenuation  effects  from  tin 


base  cloud  are  often  dominant. 
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TIME  DELAYS 


A  shock  wave  is  produced  in  the  air  by  the  explosion  of  a  muni¬ 
tion.  We  have  scaled  the  initial  radius  of  the  dust  cloud  to  the  re¬ 
gion  of  strong  shock  effects.  The  strong-shock-induced  air  flow 
fields  and  the  ejected  dust  particles  will  disrupt  the  smooth,  mean 
wind  airflow.  It  will  take  a  finite  time  before  the  mean  wind  flow 
is  reestablished  within  the  cloud  volume.  Thus  there  is  a  finite 
time  delay  before  the  dust  cloud  is  transported  horizontally  by  the 
mean  wind. 

We  define  a  characteristic  time  of 

tc  "  ^  s  ’  <2°) 

c 

where  D  is  the  characteristic  distance  and  v  is  the  characteristic 
c  c 

velocity.  For  the  base  cloud  we  take  the  characteristic  distance  to 
be  the  diameter  of  the  strong  shock  region  along  the  wind  track 
direction, 


Dc  =  2RJW  m  ,  (21) 

and  we  take  the  characteristic  velocity  to  be  the  mean  wind  velocity 
at  the  altitude  of  the  base  cloud  centroid, 

vc  =  Vw  m  s"1  .  (22) 


The  time  delay  should  be  proportional  to  the  characteristic  time.  A 
proportionality  factor  of  2  gives  the  best  fit  to  the  Dugway  transmis¬ 
sion  test  data  of  Reference  4.  The  time  delay  for  the  base  cloud  is 
then 


tDB  =  2tc 


4RTW 


W 


(23) 


The  characteristic  distance  for  the  main  cloud  is  the  same,  but 
we  take  the  characteristic  velocity  as 


where 


=  mean  wind  speed  at  a  characteristic  main  cloud  altitude 
of  10  meters  (m  s  *) 

VR  =  initial  buoyant  rise  velocity  of  the  main  cloud  (ms  1). 

The  main  cloud  time  delay  is 


=  ?t  = 
ny  “  o 


4R 

T1V 


—2  2 

V  +  V 
KM  R 


VERTICAL  WIND  PROFILE 


(25) 


The  mean  wind  speed  model  utilizes  the  following  two  assumptions. 
The  direction  of  the  wind  does  not  change  with  altitude.  The  wind 
speed  as  a  function  of  altitude  has  the  form 


V(Z)  =  V(ZR) 


(26) 


where 

V(Z)  =  m^an  wind  speed  at  altitude  Z  (m  s  ^) 

V(ZRJ  =  mean  wind  speed  at  reference  altitude  ZR  (m  s  ) 
P^,  =  power  law  exponent  of  the  vertical  wind  profile. 


The  current  model  for  horizontal  transport  of  the  dust  clouds  is 
very  simple.  After  the  tine  delays  the  lightest  particles  in  the 
clouds  are  transported  horizontally  with  a  constant  velocity.  For 
the  base  cloud,  the  constant  norizontal  velocity  is  the  mean  wind 
speed  at  the  altitude  tf  the  base  cloud  centroid.  The  main  cloud 
rises  over  a  rang  ot  altitudes.  We  choose  the  mean  wind  speed  at 
10  meters  as  the  main  cloud  horizontal  transport  velocity. 


MAIN  CLOUD  FALLOUT 

In  the  fallout  model  of  Reference  1,  we  assumed  the  altit  a.  of 
the  centroid  of  size  group  i  was  given  by 


Z  .  (t)  =  Z  (t)  -  V_,.t  m 
civ  '  cov  *  Ti 


(27) 
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where 


Zci(t)  =  altitude  of  centroid  of  size  group  i  at  time  t  (m) 

Zco(t)  =  altitude  of  centroid  of  massless  particles  (altitude 
of  ideal  buoyant  cloud)  (m) 

V„.  =  terminal  fall  velocity  under  gravity  of  size  group  i 

1  - 1 

(ms  j . 

This  formulation  implicitly  assumes  that  the  size  group  experiences 
the  full  upward  flow  of  the  buoyant  rise.  But  as  the  larger  particles 
fall  below  the  buoyant  cloud,  the  upward  flow  velocities  experienced 
by  the  particles  decrease.  Hence  the  particles  will  fall  faster  than 
our  simple  formulation. 

As  a  first-order  correction  to  allow  larger  particles  to  fall 
faster,  we  reformulate  the  fallout  to  be 

Zci(t)  ■  Zco(o)  *  [Zco">  -Zco(0>][‘  ‘)J  -VTil  * 

(28) 

where 


Z  (0)  =  initial  altitude  of  buoyant  cloud  (m) 

C0  -1 
=  initial  rise  velocity  of  buoyant  cloud  (ms  ). 

This  formulation  assumes  that  the  particles  experience  only  a  fraction 
of  the  upward  flow  velocity.  The  larger  the  particle  terminal  velocity 
is,  the  smaller  the  fraction.  The  smallest  particles  with  essentially 
zero  terminal  fall  velocity  experience  the  full  cloud  upward  flow 
velocity,  while  the  large  particles  with  terminal  fall  velocities 
equal  to  or  greater  than  the  initial  rise  velocity  do  not  experience 
any  of  the  upward  flow  velocity. 
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SECTION  3 

COMPARISONS  WITH  TEST  DATA 


INTRODUCTION 

All  of  the  models  developed  in  this  report  and  Reference  1  have 

been  incorporated  into  a  munition  dust  cloud  and  propagation  computer 

code  called  ASL-DUST.  In  this  section  we  assess  the  validity  of  the 

models  in  ASL-DUST  by  comparing  the  code  results  with  experimental 

test  data.  In  Reference  1  we  compared  the  model  results  for  the  main 

cloud  vertical  and  lateral  dimensions  as  a  function  of  time  with  the 

test  results  from  Dugway  (Reference  4)  and  Fort  Sill  (Reference  5).  In 
• 

this  report  wc  compare  the  model  size  and  altitude  predictions  with 
the  more  recent  DIRT-1  test  results  (References  6  through  11).  But 
our  primary  goal  is  to  compare  the  end  product  of  the  code  with  the 
test  results.  That  is,  we  want  to  compare  the  code  predictions  for 
electro-optical  propagation  through  the  clouds  with  the  actual  measured 
results.  This  propagation  comparison  is  the  ultimate  validity  test  for 
a  munitions  dust  code.  We  compare  code  results  with  the  infrared  trans¬ 
mission  data  from  Dugway  and  with  the  visible,  infrared,  and  milli¬ 
meter  transmission  data  from  DIRT-1.  We  do  not  have  sufficient  test 
site  data  for  Fort  Sill  to  make  code  predictions  to  compare  with  the 
transmission  data. 

MODEL  COMPARISONS  WITH  DIRT-1  TEST  DATA 
Site  Geometry 

DIRT-1  was  a  coordinated  program  to  produce  information  c  direct 
use  to  the  E-0  sensor  and  obscuration  modeling  communities.  '•  pure  2, 
taken  from  Reference  6,  shows  the  test  site  layout.  The  path  length 
from  the  north  site  to  the  south  site  is  2  kilometers.  The  test  area, 
where  explosives  were  detonated  and  artillery  shells  impacted,  is  100 


k 
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by  300  meters.  Table  2,  also  from  Reference  6,  gives  the  event  schedule 
The  explosives  for  the  events  using  three  TNT  charges  or  three  static 
projectiles  were  placed  on  a  line  perpendicular  to  the  site  center- 
line,  one  explosion  on  the  centerline  and  one  on  each  side  15  meters 
away. 

We  place  our  model  coordinate  system  with  origin  at  the  ground  sur¬ 
face  at  the  midpoint  of  the  site  centerline  with  the  Y  axis  along  the 
centerline  pointing  north,  the  X  axis  perpendicular  to  the  centerline, 
and  the  Z  axis  vertical.  The  centerline  (and  Y  axis)  lies  14.88°  west 
of  north.  Table  3  shows  the  transmitter-receiver  coordinates  assumed 
for  the  model.  i'ach  set  of  coordinates  defines  a  sight  path,  which  ex¬ 
tends  from  the  transmitter  to  the  receiver.  These  coordinates  were 
calculated  from  the  sight  path  data  given  in  Reference  6. 

Crater  Scaling  Factors 

Two  critical  parameters  required  by  any  dust  obscuration  model  are 
the  amount  and  type  of  material  lofted  into  the  air  by  a  munition  det¬ 
onation.  The  amount  of  material  can  be  related  to  the  size  of  the  ap¬ 
parent  crater  left  in  the  soil  by  the  explosion.  The  model  assumes 
that  the  scaling  relation  between  the  volume  of  the  apparent  crater 
and  the  yield  of  the  munition  is 

VA-  SAC"1,U1  11,3  »  (29) 

where 

3 

V.  =  volume  of  apparent  crater  (m  ) 

A  3-1  111 

=  apparent  crater  volume  scaling  factor  (m  (lb  TNT)  '  ) 

W  =  total  munition  yield  (lb  TNT). 

The  scaling  factor  depends  on  the  type  of  soil,  type  of  munition,  and 
depth  of  burst.  In  Reference  1  we  presented  illustrations  o  g  the 
scaling  factor  as  a  scaled  function  of  depth  of  burst  for  st-.c.  L  dif¬ 
ferent  soil  types. 

The  data  gathered  in  DIRT-1  were  the  most  complete  of  any  munition 
dust  test  series  to  that  date.  An  extensive  crater  measurement  program 


Date 

(Oct 

1978) 

Time 

(MOT) 

TABLE  2. 

Event 

EVENT 

Size 

(1b) 

SCHEDULE  FOR  DIRT-1 

Comaents 

2 

0707 

A-l 

3 

Three  1-lb  TNT  charges 

2 

0718:30 

A-2 

45 

Three  15-lb  TNT  charges 

2 

0729 

A- 3 

90 

Three  30-lb  TNT  charges 

2 

0739:50 

A-4 

90 

Three  30-lb  TNT  charges 

3 

0710 

8-1 

45 

Three  15-lb  TNT  charges 

3 

0723 

B-2 

45 

Three  15-lb  TNT  charges 

3 

0731 

B-3 

90 

Three  30-lb  TNT  charges 

3 

0741:10 

B-4 

90 

Three  30- 1 b  TNT  charges 

3 

0751  :20 

B-5 

180 

Three  60-lb  TNT  charges 

3 

0801:25 

8-6 

180 

Three  60-lb  TNT  charges 

3 

0815 

B-7 

360 

Three  120-lb  TNT  charges 

3 

0828:25 

B-8 

360 

Three  120-lb  TNT  charges 

5 

0712-50 

C-l 

2100 

One  hundred  and  forty  15-lb  charges 

6 

0717:25 

0-1 

36 

Twelve  3-lb  TNT  charges 

6 

0730:40 

0-2 

72 

Twelve  6-lb  TNT  charges 

6 

0748:55 

0-3 

144 

Twelve  12-lb  TNT  charges 

6 

0807 

0-4 

288 

Twelve  24-lb  TNT  charges 

10 

0900 

C-2 

2100 

One  hundred  and  forty  15- lb  TNT  charges 

11 

0747:40 

E-l 

45 

Three  155-mm  projectiles  (TNT),  static 

11 

0757:10 

E-2 

45 

Three  155-mm  projectiles  (TNT),  static 

11 

0821  -.45 

E-3 

180 

Twelve  155-mm  projectiles  (TNT),  static 

11 

0829:05 

E-4 

180 

Twelve  155-mm  projectiles  (TNT),  static 

12 

0737:10 

E-1 

60 

Four  155-mm  projectiles  fired  simultaneously 

12 

0748:10 

F-2 

180 

Twelve  155-mm  projectiles,  rapid  fire 

12 

0805:35 

F-3 

180 

Twelve  155-mm  projectiles,  rapid  fire 

12 

0823:10 

F-4 

120 

Eight  155-mm  projectiles,  rapid  fire 

13 

0738:10 

F-5 

120 

Eight  155-mm  projectiles,  rapid  fire 

13 

0753:35 

F-6 

120 

Eight  155-mm  projectiles,  rapid  fire 

13 

0808:20 

F-7 

120 

Eight  155-mm  projectiles,  rapid  fire 

13 

0826:15 

F-8 

45 

Three  155-mm  projectiles  fired  simultaneously 

14 

0739:50 

E-5 

45 

Three  155-mm  projectiles,  2- ft  deep,  static,  . 

14 

0756:00 

E-6 

45 

Three  155-mm  projectiles,  2- ft  deep,  static,  i 

14 

0808:45 

E-7 

45 

Three  155-mm  projectiles,  1-ft  deep,  static,  . 

14 

0823:55 

E-8 

45 

Three  155-mm  projectiles,  1-ft  deep,  static,  ! 

14 

0836:10 

E-9 

45 

Three  155-mm  projectiles,  surface,  static,  30' 

14 

0846:10 

E- 1 0 

45 

Three  155-mm  projectiles,  surface,  static,  30 

14 

1056 

G-l 

• 

Diesel  oil,  rubber,  gasoline  fire 

14 

1105:30* 

G-l 

- 

Diesel  oil,  rubber,  gasoline  fire 

14 

1108:50* 

G-l 

- 

Diesel  oil,  rubber,  gasoline  fire 

14 

1111:05* 

G-l 

- 

Diesel  oil,  rubber,  gasoline  f’re 

14 

1113:55* 

G-l 

- 

Diesel  oil,  rubber,  gasoline  fire 

14 

1116:40* 

G-l 

- 

Diesel  oil,  rubber,  gasoline  fire 

14 

1119* 

G-l 

- 

Diesel  oil,  rubber,  gasoline  fire 

14 

1122:55* 

G-l 

- 

Diesel  oil,  rubber,  gasoline  fire 

14 

1125:15* 

G-l 

• 

Diesel  oil,  rubber,  gasoline  fire 

14 

1127:15* 

G-l 

- 

Diesel  oil,  rubber,  gasoline  fire 

14 

1129:35* 

G-l 

- 

Diesel  oil,  rubber,  gasoline  fire 

14 

1132:45* 

G-l 

- 

Diesel  oil,  rubber,  gasoline  fire 

*Conmand  for  range  cameras  on;  helicopter  approximately  10  to  15  seconds  from  smoke  cloud. 


SOURCE:  REFERENCE  6 
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Instrument 

Transmitter 

Receiver 

sa 

mm 

z 

B 

D 

z 

SRI  ruby  lidar,  0.69  pm 

11.4 

-1000 

3.2 

-9.7 

1000 

4.66 

ASL  CO2  lidar,  10.6  pm 

12.3 

-1000 

3.2 

-10.7 

1000 

3.2 

NRL  fvter  transmissometer, 
0.55  m  and  10.35  urn 

1.2 

-1000 

1.7 

1.2 

1000 

1.7 

ASL  millimeter  trans- 
mi  ssometer,  94  GHz  and 

140  GHz 

-5.9 

-1000 

9.1 

1000 

2.9 

was  carried  out.  Table  4  shows  the  average  depth  and  average  diameter 
of  the  apparent  craters.  These  data  wc-e  taken  from  Reference  7.  If 
we  assume  that  the  apparent  crater  profile  is  a  section  of  a  sphere, 
then  the  volume  of  the  apparent  crater  is 


T7d  /  3D2  ,2 

a  =  T"  '  — ■.  +  d 

A  6  \  4 


3 


m 


(30) 


where 

d  =  apparent  crater  depth  (in) 

D  =  apparent  crater  diameter  (m) . 

Knowing  the  apparent  crater  volume  and  the  event  yield,  we  can  solve 
for  the  apparent  crater  volume  scaling  factor.  The  scaling  factor, 

SAC  =  ^rtn  i"3  (lb  TNT)-1-111  ,  (31) 

is  also  given  in  Table  4.  We  take  these  scaling  factor  values  as  the 
inputs  to  the  dust  model. 

Particle  Size  Distributions 

The  attenuation  due  to  dust  particles  at  a  point  along  the  sight 
path  is  determined  by  the  dust  mass  density  at  the  point,  the  complex 
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TABLE  4,  APPARENT  CRATER  DIMENSIONS  AND  SCALING  FACTORS 


Event 

Yield 

W  (lb  TNT) 

Depth 

d(m) 

Diameter 

D(m) 

Scaling  Factor 

SAC  (m3[lb  TNT]"1*111) 

A1 

1 

ssi 

0.85 

0.0869 

A2 

15 

1.75 

A3 

30 

K£j 

2.23 

0.0322 

A4 

30 

0.60 

2.23 

0.0293 

B1 

15 

0.49 

1.43 

0.0225 

B2 

15 

0.54 

1.40 

0.0246 

B3 

30 

0.59 

1.98 

0.0233 

B4 

30 

0.62 

1.85 

0.0219 

B5 

60 

0.86 

3.02 

0.0361 

B6 

60 

0.95 

3.88 

0.0642 

B7 

120 

1.02 

3.67 

0.0291 

B8 

120 

0.96 

3.87 

0.0299 

Cl 

15 

0.57 

1.66 

0.0352 

C2 

15 

0.55 

1.76 

0.0373 

01 

3 

0.38 

1.03 

0.0552 

02 

6 

0.43 

1.15 

0.0362 

D3 

12 

0.52 

1.49 

0.0333 

04 

24 

0.57 

2.00 

0.0290 

El 

«15  (155  mm 
shell ) 

0.49 

1.52 

0.0250 

E2 

0.49 

1.40 

0.0117 

E3 

■  0.52 

1.40 

0.0234  ( 

E4 

:  0.52 

1.44 

0.0245 

E5 

0.74 

2.77 

0.120  ! 

E6 

0.75 

2.59 

0.109 

E7 

0.66 

2.36 

0.0785 

E8 

0.66 

2.26 

0.0726 

E9 

0.50 

1.79 

0.0343 

E10 

0.34 

1.65 

0.0190 

Fa 

0.37 

1.44 

0.0162 

aAverage  values  for  events  F1-F8. 
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index  of  refraction  of  the  dust  particles,  and  the  dust  particle  size 
distribution.  Due  to  dust  particle  fallout  and  diffusion,  the  particle 
size  distribution  varies  as  a  function  of  position  within  the  cloud  and 
time  after  burst.  By  following  the  various  size  groups  as  a  function  of 
time,  the  model  can  calculate  the  particle  size  distribution  at  a  point 
provided  the  size  distribution  of  the  initial  particles  lofted  is  known. 
At  present  there  is  .10  adequate  method  of  predicting  the  initial  lofted 
particle  size  distribution  from  the  standird  physical  soil  measurements. 
More  munition  tests  n  different  soil  ■  ypes  are  needed  before  a  general 
size  distribution  model  can  be  formulated  with  any  confidence. 

One  standard  soil  procedure  is  to  use  sieve  and  hydrometer  analysis 
on  bulk  soil  samples  to  determine  the  soil  grain  distribution.  In  this 
analysis  all  of  the  grains  are  separated  from  each  other  before  measure¬ 
ment.  Figure  3,  reproduced  from  Reference  7,  shows  a  cumulative  grain 
mass  distribution  for  the  DIRT-1  soil.  Figure  4  shows  this  data  replot¬ 
ted  on  log-log  paper.  Also  shown  on  this  figure  arc  the  cumulative 
masses  calculated  from  three  standard  number  size  probability  distribu¬ 
tions,  two  power  laws  and  one  log-normal.  We  see  that  a  reasonable  fit 
to  the  soil  grain  Sxze  distribution  would  be  a  oower  law  distribution 
with  exponent  p  =  5.05  for  grains  up  to  a  diameter  of  60  microns  joined 
to  a  log-normal  distribution  (a^  =7.1  microns,  '>  -  ?.65)  for  particles 
greater  than  60  microns. 

Can  wo  use  the  measured  soil  grain  size  distribution  as  the  model 
initial  lofted  dust  particle  distribution?  We  can  if  the  explosion 
dynamics  simply  separate  and  loft  the  in  situ  soil  grains.  However, 
if  in  the  cratering  and  lofting  processes  some  grains  are  fractured 
and/or  remain  stuck  together,  then  the  lofted  initial  size  distribution 
will  differ  from  the  soil  grain  distribution. 

An  experiment  to  measure  the  actual  particle  size  distril  tions 
within  the  explosively  produced  dust  clouds  was  fielded  in  D 1 1<  .. 

This  experiment  consisted  of  three  particulate  spectrometers  moiu.t'  J 
on  an  airborne  platform;  the  platform  was  suspended  beneath  a  helicop¬ 
ter,  which  maneuvered  it  through  the  dust  clouds  at  various  altitudes 
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Figure  4.  Cumulative  grain  mass  distribution  versus  grain  size. 


and  times  after  burst.  Figure  5  (from  Reference  6)  shows  the  particle 
size  distribution  measured  at  10.6  meters  above  ground  level  about  2 
minutes  after  completion  of  an  eight-round  barrage  of  155-mm  super¬ 
quick  point  detonating  projectiles  (event  F6) .  Note  that  this  figure 
shows  the  number  size  distribution,  in  contrast  to  Figures  3  and  4 
which  show  mass  cumulative  probabilities. 

On  the  log- log  graph  of  Figure  5  a  power  law  distribution  plots  as 
a  straight  line.  Fitting  a  straight  line  through  the  data  points  gives 
a  power  law  exponent  of.  about  3.4.  This  value  is  intermediate  between 
the  power  law  fits  to  the  soil  grain  data  of  exponents  3.65  for  grains 
less  than  40  microns  in  diameter  and  3.2  for  grains  bigger  than  40 
microns . 

Integrating  the  size  distribution  and  normalizing,  we  have  the 
cumulative  probability  size  distribution  shown  in  Figure  6  plotted  on 
lognormal  probability  paper.  The  cumulative  distribution  from  any 
lognormal  size  probability  distribution  plots  on  this  probability  paper 
as  a  straight  line.  The  fit  to  the  data  shown  in  Figure  6  is  a 
lognormal  size  probability  distribution  with 

a  =  3 . 1  pm 
m 

S  =  2.34 

The  airborne  particle  data  was  taken  about  2  minutes  after  burst. 
Fallout  would  have  reduced  the  airborne  populations  of  dust  particles 
larger  than  about  50  to  100  microns.  But  airborne  particles  less  than 
about  10  microns  should  have  been  virtually  unaffected  by  fallout  at 
this  time.  The  airborne  data  show  relatively  less  small  particles 
than  the  soil  grain  data.  Assuming  both  sets  of  data  are  accurate, 
the  difference  suggests  that  not  ail  of  the  smaller  soil  grains  are 
being  separated  during  the  cratering  and  lofting  dynamics. 

Because  of  the  fallout,  the  airborne  data  are  of  limited  use  in 
predicting  the  initial  size  distribution  of  particles  greater  than 
about  50  to  100  microns  in  diameter.  Using  the  code,  we  input  a  number 
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Figure  5.  Airborne  particle  size  distribution. 
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Figure  6.  Airborne  cumulative  size  distribution. 

of  size  distributions  for  the  larger  particles  and  compared  the  pre¬ 
dicted  transmissions  with  test  data.  The  log-normal  fit  to  the  larger 
particles  in  the  grain  size  data,  in  Figure  4,  produces  a  reasonable 
fit  to  the  data.  For  the  model  we  assume  a  bimodal  dust  distribution. 
We  take  mode  A,  the  small  particles,  to  have  the  log-normal  fit  to  the 
airborne  data,  and  take  mode  B,  the  larger  particles,  to  have  the  log¬ 
normal  fit  to  the  soil  grain  data.  Thus 


Mode  A :  a  =  3 . 1  ym 
m 

S  =  2.34 

Mode  B:  a  *  7.1  um 
m 

S  =  2.65 

Meterological  Conditions 

The  wind  is  often  the  most  important  parameter  in  determining  the 
strength  an  duration  of  degradation  of  the  electro-optical  signals  in 
a  dusty  environment.  The  wind  can  blow  a  dust  cloud  into  or  out  of  a 
sensor  sight  path.  The  stronger  the  wind,  the  faster  the  atmospheric 
diffusion  will  disperse  the  dust  cloud. 

In  DIRT-1  two  meteorological  data  collection  sites  were  employed, 
one  at  each  end  of  the  2  km  optical  path  (see  Figure  2).  At  the  south 
site  the  wind  direction  and  speed  were  measured  at  1,  2,  and  4  meters 

above  ground  level.  At  the  north  site  only  the  2  meter  winds  were 

recorded.  The  wind  data  were  averaged  for  both  a  5-minute  and  a  10- 
min"te  period  starting  at  the  beg;nning  of  each  test  event.  The  winds 
during  all  test  events  were  light  and  extremely  variable.  lhe  wind 

data  given  in  Reference  6  show  that  the  average  wind  direction  at  one 

site  sometimes  changed  direction  by  more  than  45  degrees  from  one  10- 
minute  averaging  period  to  the  nc.  l .  Moreover,  the  wind  direction  at 
the  two  ends  of  the  sight  path  at  the  same  time  were  almost  uncorrelated, 
often  being  180  degrees  out  of  phase. 

For  the  dust  model  we  need  to  specify  the  mean  wind  speed  at  a 
reference  altitude,  the  exponent  of  the  vertical  profile,  and  the  wind 
direction.  We  take  2  meters  as  the  reference  altitude,  and  use  the  10- 
minute  averaged  speed  at  2  meters  at  the  south  site  as  the  mean  wind 
speed.  The  wind  direction  at  the  midpoint  of  the  sight  path  is  diffi¬ 
cult  to  estimate  because  of  the  extreme  variability.  We  make  a  f irst 
estimate  by  averaging  the  north  and  south  site  directions.  Wc  nen 
modify  the  direction  to  produce  the  best  agreement  between  model  and 
test  transmission  results.  For  each  event  we  calculate  an  aveiage 
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vertical  profile  exponent  using  the  south  site  1-,  2-,  and  4-meter 
data.  We  then  average  all  exponents  to  come  up  with  one  average  ex¬ 
ponent  of  0.19,  which  we  use  for  &11  events.  Table  5  shows  the  mean 
wind  speeds  and  directions  chosen  for  the  dust  model  inputs. 

Other  Model  Parameters 

From  the  soil  measurements  of  Reference  7,  we  take  the  average 
moisture  fraction  of  the  soil  to  be  8.4  percent  and  the  average  in  situ 

_3 

soil  density  to  be  1.6  g  cm  .  Table  6  shows  the  indices  of  refrac¬ 
tion  used  in  the  model.  The  visible  and  infrared  values  are  the  values 
recommended  in  Reference  6,  the  carbon  values  are  from  Reference  12,  and 
the  millimeter  values  were  chosen  to  provide  agreement  between  model  and 
test  transmission  results. 

At  its  present  stage  of  development,  ASL-DUST  does  not  include 
multiple-burst  interactions.  Although  the  model  allows  up  to  10 
bursts,  the  size,  position,  and  internal  mass  distribution  of  each 
dust  cloud  are  calculated  independently  of  all  other  clouds.  The 
total  propagation  effect  is  simply  taken  as  the  sum  of  the  individual 
effects.  If  the  dust  clouds  from  the  multiple  bursts  are  separated 
sufficiently  so  there  is  little  mutual  overlap,  noninteraction  is  a 
good  assumption.  However,  if  the  multiple  bursts  are  spaced  close 
enough  in  time  and  space,  then  the  rising  dust  clouds  will  intersect; 
the  air  flow  fields  will  interact  and  rearrange  to  form  one  larger 
merged  cloud  from  the  individual  clouds.  The  single  merged  cloud  will 
have  different  transport,  diffusion,  and  propagation  properties  than 
the  sum  of  the  individual  clouds. 

All  of  the  DIRT-1  test  events  were  multiple  bursts.  For  all  but 
the  smallest  yields,  the  individual  clouds  rapidly  merged  to  form  a 
single  cloud.  For  the  smallest  yields  we  input  the  individual  bursts. 
For  the  other  yields  we  simulate  the  single  merged  cloud  in  the  model 
by  the  following  method.  We  input  only  a  single  charge  at  the  center 
of  the  multi-charge  layout;  the  yield  of  the  single  charge  is  the  sum 
of  the  individual  yields.  When  the  individual  clouds  first  merge,  the 


TABLE  6.  COMPLEX  INDICES  OF  REFRACTION 


Wavelength  or 
Frequency 

Material 

Dust 

Carbon 

0.55  ym 

1.525  -  i  0.005 

1.53  -  i  0.005 

10.6  ym 

1.65  -  i  0.14 

1.65  -  i  0.14 

94  GHz 

1.45  -  i  0.033 

2  -  i  1 

merged  single  cloud  is  larger  than  the  cloud  from  an  equivalent  single 
burst.  Also  the  rise  rate  of  the  merged  cloud  is  smaller  than  the 
equivalent  single-burst  cloud.  We  roughly  account  for  both  of  these 
effects  in  the  model  by  increasing  the  entrainment  factor,  a,  from 
its  nominal  value  of  1  to  a  value  of  2.  This  increased  a  factor  will 
cause  the  equivalent  single  cloud  to  entrain  ambient  air  more  rapidly. 
The  increased  entrainment  will  cause  the  cloud  to  grow  more  rapidly  in 
size  and  to  rise  at  a  slower  rate,  thereby  more  closely  simulating  the 
behavior  of  the  merged  cloud. 

In  general,  the  other  model  parameters  not  mentioned  were  set  at 
the  model  default  values. 

Size  and  Altitude  Comparisons 

When  Reference  1  was  prepared,  only  the  Fort  Sill  and  Dugway  test 
data  were  available;  we  did  not  have  any  cloud  altitude  test  data  to 
compare  against  the  model  predictions.  The  DIRT-1  test  data  does  in¬ 
clude  cloud  altitude;  Figures  7  through  13*  show  the  ASL-DUST  model 
comparisons  with  the  test  data  for  the  B-series  of  events.  The  test 
data  shown  are  for  the  centroid  of  the  whole  main  cloud.  The  model 
data  shown  are  for  the  centroid  of  the  size  group  of  the  smallest 
particles.  Until  fallout  becomes  significant  at  later  times,  these 
two  centroids  should  be  approximately  equal.  At  later  times  we  expect 
the  whole  cloud  centroid  to  fall  below  the  smallest  particles  centroid. 

*Figures  7  through  55  are  on  pages  45  through  84. 


In  general  the  intermediate  and  later  time  comparisons  are  as  expected; 
the  two  centroids  are  approximately  equal,  and  then  the  test  data  cen¬ 
troid  falls  below  the  model  centroid.  In  the  first  few  seconds  the 
test  data  centroid  is  higher  than  the  model  centroid.  The  model  alti¬ 
tude  is  that  of  a  buoyantly  rising  cloud.  However,  for  about  the  first 
second,  the  cloud  rise  is  controlled  by  the  shock-induced  air  flow 
fields.  These  predominantly  upward  flowing  fields  carry  the  cloud  up¬ 
ward  one  or  two  cloud  diameters.  The  shock-induced  rise  forces  die 
away  and  the  buoyant  r^se  fo^ce  controls  the  rise  history  thereafter. 

This  first-second  nonbuoyant  rise  regime  should  be  scalable  from  the 
shock  and  initial  cloud  dimensions.  A  model  for  this  effect  will  be 
abided  to  ASL-DUST  at  a  later  date. 

Figure  14  shows  the  cloud  geometry  of  event  B7  with  the  model  pre¬ 
dictions  superimposed.  This  event  is  three  120-lb  charges  in  a  line 
15  meters  apart.  The  early  time  cloud  is  thus  wider  than  the  cloud  from  a 
single  560- lb  charge.  In  the  model  we  simulate  this  three-charge  event 
with  a  single  360-lb  charge  with  an  enhanced  entrainment  coefficient. 

The  geometry  comparisons  show  the  expected  earlv  time  difference,  but 
the  comparisons  are  quite  good  from  about  10  seconds  on.  The  compari¬ 
sons  at  20  seconds  illustrate  the  differences  in  the  two  centroids. 

The  model  centroid  and  geometry  are  a  good  fit  to  t.ie  upper  cloud, 
where  wo  expect  the  smaller  particles  to  reside.  In  the  model,  the 
lower  part  of  the  cloud  is  accounted  for  by  the  fallout  of  the  larger 
size  groups. 

Attenuation  Comparisons 

The  DIRT-1  B-series  of  events  is  particularly  valuable  for  model 
comparisons.  The  events  have  a  fairly  large  yield  range,  from  three 
15-lb  charges  to  three  120- lb  charges.  Test  data  were  taken  on  crater 
size,  cloud  size  and  altitude,  and  cloud  attenuation  at  visible,  infra¬ 
red,  and  millimeter  wavelengths.  All  these  test  data  provide  a  good 
test  for  any  dust  cloud  model.  The  wide  range  of  wavelengths  in  the 
attenuation  measurements  provide  both  early  and  late  time  comparisons 
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as  well  as  the  spectral  sensitivity  comparisons.  The  significant  milli¬ 
meter  attenuation  generally  lasts  less  than  about  10  seconds;  the  vis¬ 
ible  and  infrared  attenuation  can  last  for  hundreds  of  seconds.  Thus 
the  millimeter  data  provide  a  data  base  for  validation  of  the  early 
time  model,  while  the  visible  and  infrared  data  provide  a  late  time 
data  base. 

Figures  15  through  25  show  the  data  (Reference  8)  and  model  compar¬ 
isons  at  94  GHz  for  events  B4  through  B8,  and  events  E5  through  E10 
For  the  B  series  the  chart  recorder  time  constant  was  400  ms,  and  for 
the  E  series  it  was  changed  to  40  ms.  For  the  B  series,  the  maximum 
magnitude  of  the  model  attenuation  predictions  and  the  time  history 
are  quite  similar  to  the  measured  test  data.  Because  of  the  model  as¬ 
sumption  that  the  initial  cloud  *orms  instantaneously  at  burst  time, 
the  model  attenuation  peaks  at  burst  time  and  declines  thereafter. 

The  real  cloud  takes  a  finite  time  to  form  and  build  up  to  maximum  at¬ 
tenuation.  The  figures  show  that  the  test  data  attenuation  reaches  a 
maximum  about  1.5  seconds  after  burst  time  for  the  B  series  and  slightly 
less  than  a  second  for  the  E  series.  Part  of  the  delay  is  artificial 
and  is  due  to  the  chart  recorder  time  constant,  and  part  is  real  and 
is  due  to  the  cloud  formation  physics.  The  real  part  of  the  delay  is 
probably  about  one-half  to  three-quarters  of  a  second.  We  could  im¬ 
prove  the  agreement  between  model  and  data  in  the  first  second  after 
burst  by  replacing  the  model  instantaneous  cloud  formation  with  a  finite 
formation  tir.e.  At  present  this  level  of  detail  is  not  considered 
necessary.  If  at  a  later  time  som^  E-0  system  is  sensitive  to  these 
first-second  details,  then  the  model  improvement  can  be  implemented. 

In  the  E  series  of  tests,  note  the  prominent  spikes  in  the  test 
data.  The  shorter  time  constant  of  the  E  series  allows  these  scintilla¬ 
tions  to  be  resolved.  Dust  clouds  are  turbulent  and  the  spatial  dis¬ 
tribution  of  the  lofted  dust  is  inhomogeneous.  The  model  predictions 
are  for  the  mean  attenuation;  there  the  agreement  is  quite  reasonable. 

The  E  series  of  tests  is  a  good  illustration  of  the  importance  of 
the  initial  cloud  dust  mass  loading  and  also  verifies  the  ASL-DUST 
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model  of  this  important  parameter.  In  our  model  we  assume  the  initial 
lofted  dust  mass  is  proportional  to  the  soil  mass  that  was  in  the  ap¬ 
parent  crater  volume  left  by  the  explosion: 


where 


‘CM 

PG 

FH20 


=  total  mass  of  dust  grains  lofted  (g) 

=  fraction  of  mass  in  apparent  crater  that  is  lofted 

-3 

-  bulk  density  of  tr.  sztu  soil  (g  cm  ) 

=  soil  moisture  fraction  by  weight 
=  volume  of  apparent  crater  (cm'*) . 


The  only  model  input  parameters  that  change  from  event  to  event  in 
the  E  series  are  the  apparent  crater  volume  and  the  wind  speed  and 
direction.  But  the  winds  are  so  light  that  they  have  essentially  no 
effect  upon  the  calculated  attenuation  during  the  few  seconds  of  sig¬ 
nificant  attenuation.  Hence  all  model  differences  from  event  to  event 
are  due  to  the  changing  crater  volumes.  These  events  were  all 
static  detonations  of  three  lS5-mm  projectiles,  tilted  at  30  degrees. 
In  events  E5  and  F6  the  projectiles  were  buried  2  feet  deep;  in  F.7  and 
E8  they  were  1  foot  deep;  and  in  E9  and  E10  they  were  at  the  surface. 
As  expected,  the  apparent  crater  volume  increased  with  burial  depth. 

Assuming  the  apparent  crater  was  a  segment  of  a  sphere,  the  average 

3  3 

crater  volumes  were  2.3  m  for  E5  and  E6,  l.S  m  for  E7  and  E8,  and 

3 

0.54  m  for  E9  and  E10.  Both  the  calculated  and  observed  attenuations 
decrease  with  decreasing  crater  volumes. 


Figures  26  through  35  show  the  data  (Reference  11)  and  the  model 
predictions  for  the  NRL  measured  transmissions  at  visible  (0.55  pm) 
and  infrared  (10.35  pm)  wavelengths  for  events  B1  through  B8  and  El 
and  E2.  An  interesting  phenomenon  is  exhibited  by  the  data  and  mir¬ 
rored  by  the  model  predictions.  In  passing  through  particulate 
regions,  normally  the  longer  wavelength  signal  experiences  less  at¬ 
tenuation  (more  transmission)  than  the  3horter  wavelength  signal.  But 
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both  model  and  data  show  an  inversion  of  this  normal  relationship  on 
the  recovering  signal;  the  infrared  is  attenuated  mort  than  the  visi¬ 
ble.  This  rather  unexpected  result  is  due  to  the  dust  particle  size 
distribution. 

Figure  36  shows  the  Mie  extinction  efficiencies  for  the  visible 
and  infrared  wavelengths  as  a  function  of  dust  particle  diameter. 

For  particles  less  than  about  5  microns  in  diameter,  the  visible  wave¬ 
length  has  a  higher  extinction  than  the  infrared.  But  for  particles 
in  the  size  range  from  5  to  about  50  microns  (in  the  region  of  the 
infrared  extinction  peak),  the  infrared  has  a  higher  extinction.  For 
larger  particles,  both  wavelengths  approach  the  large  particle  extinc¬ 
tion  limit  of  2;  the  extinctions  are  essentially  equal. 

Thus  the  relative  visible/infrared  extinction  depends  on  the  par¬ 
ticulate  size  distribution.  The  carbon  particles  in  the  dust  cloud 
are  modeled  as  a  log-normal  size  distribution  with  a  mean  diameter  of 
0.5  microns  and  standard  deviation  parameter  of  2.  For  these  small 
particles,  the  visible  is  considerably  more  attenuated  than  the  infra¬ 
red.  And  indeed  the  attenuations  measured  in  the  oil  fire  event  (where 
we  expect  carbon  to  be  a  major  particle  constituent)  confirmed  the 
predicted  attenuation  behavior.  But  for  dust  clouds,  carbon  is  normally 
a  minor  constituent.  The  spectral  dependence  of  the  dust  clouds  is  a 
function  of  time.  At  very  early  times  when  many  large  particles  are 
still  aloft,  we  would  expect  the  vlslble/infrared  attenuations  to  be 
approximately  equal.  As  fallout  diminishes  the  larger  particles,  the 
size  distribution  of  the  particles  aloft  shifts  toward  the  smaller 
particles.  The  dust  particle  size  distributions  measured  and  assumed 
for  the  model  have  a  significant  fraction  of  particles  in  the  5  to  50  pm 
inversion  region.  Thus  after  a  short  fallout  time,  inversion  occurs; 
the  Infrared  is  attenuated  more  by  the  dust  particles  aloft  than  the 
visible.  If  we  were  to  wait  long  enough,  fallout  would  shift  the 
particle  size  distribution  far  enough  so  that  another  inversion  would 
occur  and  the  visible  would  uc  attenuated  more. 
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The  spectral  dependence  of  the  dust  cloud  is  not  only  a  function  of 
time,  but  is  also  a  function  of  position  within  the  cloud.  Smaller 
particles  tend  to  rise  to  the  top  of  the  cloud  and  the  larger  particles 
sink  to  the  bottom  of  the  cloud.  Thus  while  a  sight  path  through  the 
top  of  the  cloud  might  show  the  normal  visible/infrared  dependence,  a 
sight  path  near  the  bottom  could  simultaneously  exhibit  the  inverted 
dependence.  The  advantage  in  using  the  ASL-DUST  method  of  many  size 
groups  and  an  exact  Mie  calculation  is  that  the  model  automatically 
calculr.t  s  the  snect^al  dependence  both  as  a  function  of  time  and  posi¬ 
tion  within  the  cloud.  Simpler  models  using  an  average  mass  extinc¬ 
tion  coei^icient  cannot  model  these  spatial  and  temporal  dependencies. 

In  these  visible/infrared  model  transmission  comparisons,  the 
wind  plays  a  major  role,  in  contrast  to  the  comparisons  at  millimeter 
wavelengths.  The  recovery  time  of  the  signal  is  determined  primarily 
by  the  wind  blowing  the  cloud  out  of  the  sight  path.  The  speed  of  the 
recovery  is  determined  partly  by  the  wind  and  partly  by  the  model  mass 
distribution  and  uiffusion.  The  visible/infrarcd  spectral  dependence 
is  determined  by  the  initial  assumed  particle  sire  distribution  and 
fallout  mechanisms.  In  general  the  comparisons  arc  quite  good.  A 
real  dust  cloud  is  a  lumpy,  bumpy  thing  and  the  real  wind  is  not  steady 
in  direction.  In  events  B5  and  B7  the  signal  essentially  recovers  and 
then  another  piece  of  dust  cloud  is  blown  by  the  sight  path.  In  event 
B5,  we  chose  the  model  wind  direction  so  the  model  transmission  was 
intermediate  between  the  two  cloud  pieces.  In  B7  we  chose  to  match 
the  primary  recovery.  In  event  B8,  the  test  data  show  that  the  B8  dust 
cloud  was  initially  blown  away  from  the  line  of  sight,  and  about  8 
minutes  later  the  wind  blew  the  B8  cl  ^d  or  another  cloud  back  through 
the  line  of  sight. 

MODEL  COMPARISONS  WITH  DUGWAY  TEST  DATA 

Tn  November  and  December  of  1977  six  dust  cloud  tests  were  ,r.dueted 
at  the  Dugway  Proving  Ground.  Dust  clouds  were  generated  using  static 
TNT  charges  to  simulate  the  detonation  of  high-explosive  munitions 


rounds.  Table  7  shows  the  trial  number,  the  munition  simulated,  and 
the  amount  of  TNT  used.  In  the  previous  report.  Reference  1,  we  com¬ 
pared  the  initial  dust  cloud  size  and  the  time  histories  of  the  cloud 
width  and  hoight  with  the  model  predictions.  Now  we  will  compare  the 
measured  transmission  at  3.4  pm  with  the  model  predictions. 

TABLE  7.  LIST  OF  EXPLOSIVES  AND  SIMULATED 
ROUNDS  FOR  DUGWAY  TESTS 


Trial  Number 

Simulated  Round 

Weight  of  TNT  (lb) 

DPG-002-D1 

81-mm  mortar 

3.5 

DPG-002-D2 

4.2-inch  mortar 

5.0 

DPG-002-D3 

105-mm  tank  HEP 

6.75 

DPG-002-D4 

120-mm  mortar 

8.0 

DPG-002-D5 

155-mm  artillery  projectile 

15.0 

DPG-002-D6 

8-inch  artillery  projectile 

35.0 

In  the  DIRT-1  test  series,  the  static  filings  were  centered  about 
the  line  of  sight.  In  the  Dugway  tests,  the  dust  cloud  was  generated 
to  one  side  of  the  lines  of  sight  and  the  prevailing  wind  blew  the 
cloud  through  the  sight  paths.  Figure  37  shows  the  three  sight  paths 
and  detonation  points.  D6  is  an  the  opposite  side  of  the  site  because 
of  a  shift  in  the  wind  direction.  The  DIRT-1  comparisons  serve  as  a 
test  of  the  model  predictions  for  the  cloud  behavior  near  the  forma¬ 
tion  point.  The  Dugway  series  will  test  the  model  predictions  for  the 
cloud  parameters  after  wind  transport  over  the  ground. 

The  dust  particle  size  distributions  and  the  crater  volumes  were 
not  measured  in  the  Dugway  tests.  We  will  assume  the  Dugway  soil  to 
be  similar  to  the  DIRT-1  soil  and  use  the  DIRT-1  size  distributions 
and  crater  scaling  factors.  We  use  the  Dugway  measured  meteorology 
parameters  of  wind  directions,  wind  profile  power  exponent,  and 
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Pasquill  stability  factor.  We  choose  an  average  wind  speed  which  cor¬ 
responds  with  the  arrival  of  the  dust  cloud  at  the  most  distant  sight 
path.  Table  8  shows  the  model  input  parameters. 


TABLE  8.  MODEL  INPUT  PARAMETERS  FOR  DUGWAY  TESTS 


Input  Parameter 

Trial  Number  l 

D1 

D2 

D3 

D4 

D5 

D6 

Yield  (lb  TNT) 

3.5 

5.0 

6.75 

8.0 

15 

35 

Crater  volume  scaling 
factor  '"i3  [lb  TNT]"1,111) 

0.064 

0.06 

0.06 

0.04 

0.04 

0.032 

Wind  speed  at  8  meters 
altitude  (m  s"1) 

2.2 

4.0 

5.0 

H 

6.5 

4.7 

Wind  azimuth  (deg) 

208 

204 

171 

175 

157 

344  ; 

Wind  vertical  orofile 
power  exponent 

0.06 

0.08 

0.09 

0.07 

0.15 

0.16  , 

Pasquill  stability  factor 

C 

D 

D 

D 

D  | 

D 

Note:  Model  azimuth  =  true  azimuth  -  145°. 


Figures  38  through  55  show  the  data  (Reference  4)  and  the  model 
comparisons.  Any  individual  dust  cloud  has  a  very  irregular  structure. 
The  model  cloud  has  a  very  regular  structure  which  supposedly  repre¬ 
sents  the  average  of  a  large  number  of  real  clouds.  So  our  model  pre¬ 
dictions  can  be  expected  to  vary  about  the  measured  data  but  to  agree 
on  the  average.  For  D1  the  model  cloud  is  early  in  reaching  row  0  but 
is  late  in  reaching  row  M;  for  both  rows  the  width  of  the  cloi  J  and 
the  minimum  transmission  agree  approximately  with  the  data.  Th(  other 
tests  show  variations  between  model  and  data,  but  on  the  average  the 
comparisons  are  quite  reasonable. 
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figure  13.  Dust  cloud  centroid  altitude  for  event  B8 
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Figure  14.  Comparison  of  cloud  geometry  for  event  B7 
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Figure  18.  Attenuation  at  94  GHz  for  event  B7. 


Figure  19.  Attenuation  at  94  GHz  for  event  B8. 
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Figure  21.  Attenuation  at  94  GHz  for  event  E6. 
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Figure  22.  Attenuation  at  94  GHz  for  event  E7. 
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Figure  23. 
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Attenuation  at  94  GHz  for  event  E8. 
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Figure  25.  Attenuation  at  94  GHz  for  event  E10. 
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Figure  26.  Transmittance  at  10.35  ym  and  0.55  ym  for  event  Bl. 
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Figure  27. 


SOURCE:  REFERENCE  11 
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Transmittance  at  10.35  ym  and  0.55  ym  for  event  B2. 
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Figure  29.  Transmittance  at  10.35  um  and  0.55  ym  for  event  B4 
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Figure  30.  Transmittance  at  10.35  ym  and  0.55  ym  for  event  B5. 


Figure  31, 
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Transmittance  at  10.35  ym  and  0.55  ym  for  event  B6. 
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Figure  34. 
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Transmittance  at  10.35  urn  and  0.55  pm  for  event  El. 
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Figure  51.  Transmission  at  3.4  ym  for  event  D5„  row  0. 
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